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I. INTRODUCTION 

Before presenting the details, let me first summarize the principal features of the so-called t-J model. This model 
represents a very nontrivial model of strongly correlated fermions because of the following principal reasons: 

1. It represents a system of strongly inter-correlated itinerant fermions which transform into an antiferromagnetic 
state of localized particles (the Mott-Huhhard insulator); the classic situation takes place at the concentration 
of one fermion per single-band state (at the half filling, n = \). The half- filled state is an antiferromagnetic 
insulator modeled by the Anderson kinetic exchange. 

2. The itinerant state is represented by particles called the correlated holes or simply the holes in the Mott-Hubbard 
insulator, which do not have the ordinary fermion properties (their creation and annihilation operators do not 
obey the fermionic anticommutation rules). In other words, they cannot be represented by Landau quasiparticles 
in an exact manner, since they do not represent almost-filled band states. 

3. Magnetic interaction between the correlated itinerant particles, again the kinetic exchange, is regarded also 
as the source of real-space pairing as it is represented as taking place between the nearest neighbors. Hence, 
antiferromagnetism and this new type of paired state must be regarded on equal footing with the paired state 
(resonance- valence bond state or superconducting state). This is particularly relevant in the context of high- 
temperature superconductivity near the band filling, i.e. close to the Mott transition. 

4. Because of the nontrivial character of the kinetic- (or residual-band-) energy {the projected hopping part), it 
can become comparable or even smaller than the kinetic-exchange-energy part. In effect, magnetic polaron or 
phase-separated states or new type of spin-paired states can be formed, concomitantly with the transition to a 
localized state for small carrier concentration of holes. Thus the transition in real systems such as La2-a;Sr2,Cu04, 
YBa2Cu3 06-i-a;, Lai_2:Sra;Ti03, takes place even for a non-half-filled band configuration. Although, here the 
role of atomic disorder is probably also very important, if not crucial. 

The four mentioned above features originate from the circumstance, that the kinetic (band) energy of the fermionic 
particles in the correlated systems is relatively small and easily comparable (if not smaller) to the Coulomb interaction 
energy. Furthermore, both energies are then effectively counted on the Kelvin, rather than on the electronovolt scale 
and represent competing compensating each other, dynamical contributions to the electron states. Hence, the system 
is very susceptible to the perturbations such as thermal or atomic disorder, or the electron-lattice coupling. This 
sensitivity of the strongly correlated metallic (or magnetic- insulating) state leads to an instability with respect to 
the perturbations on the thermal scale and complicates enormously a reliable solution of the model, in addition to 
the specific situation that, unlike in the standard quantum mechanics, there is no small parameter in the model, since 
the potential (interaction) energy is at least comparable to the single-particle (band) energy. 

The structure of the paper is as follows. In Section 2 we define the limit of strong correlation. In Section 3 we 
discuss the original derivation and the meaning of t-J model. In Section 4 we rewrite the exchange part in terms of 
pairing operators and discuss this new aspect briefly, whereas in Section 5 we introduce hybrid pairing. 
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II. THE HUBBARD MODEL 



The Hubbard model (1963) [1], was devised first to describe the single-band magnetism, particularly to understand 
the so-called itinerant (Stoner-Wolfharth) ferromagnetism and criterion of its appearance in a microscopic manner. 
It is represented by the Hamiltonian, which in the second- quantization representation has the form 

H = J2' "^'^ "J'^ + C/ ^ "a ' (1) 

ija i 

where the primed summation means that we take only i ^ j terms in the hopping part. The ingenuity of this 
expression derives from the fact, that we start from the atomic (Wannier) representation rather than from the popular 
then electron-gas (band) representation of the many-electron state in solid. This representation (and actually, the 
Hamiltonian) , has been invented earlier [2] , although the meaning of the model was in the latter case quite different 
namely, to derive the antiferromagnetic superexchange interaction as an effective d—d exchange in the Mott insulators 
(such as NiO, MnO or CoO), without involving, explicitly the filled p shells of oxygen or other chalcogenide ions. 

Obviously, one should think that the Bloch {$k(i")} and the Wannier {Wi(r)} single-particle bases are physically 
equivalent, as they are related by (unitary) Fourier transform on the lattice 



<J^k(r) = W.(r)e^k.R. 



N ■ - (2) 



Likewise, the same unitary equivalence holds between the annihilation operators in the Bloch {okcr} and the Wannier 
{flio-} namely 



- VTV «ko 



(3) 



The main (and fundamental) aspect of this problem is that, depending on the value of [/, the equivalence between 
the two representations can be broken and, as U increases, the atomic representation of electronic states becomes more 
appropriate than the Bloch one. This fact expresses the so-called Mott-Huhhard or metal-insulator transition, which 
in this case, is associated with the breakdown of the global C/(l) symmetry, as in the Mott-Hubbard insulating phase 
linear momentum is in an obvious manner not a proper quantum number. Parenthetically, one should say that the 
Hubbard (1964) method of defining the localization [3] is similar to that of Mott. In essence, the Hubbard approach 
takes into account both the existence of crystal lattice, as well as the narrow-band structure of the correlated electrons. 
Explicitly, taking expectation values of the two terms in ([1]) one sees, that in the Hartree-Fock approximation and for 
one electron per atom we have the ground state energy per site in the form 



- -\Y. I II (1 - ^^-) + Un^n^ = -— + - , (4) 

where W is the bandwidth and z number of nearest neighbors {W = 2z|i|). In the paramagnetic state {ria = n-a = 1/2) 
therefore, the two energies compensate each others when U = W i.e. when band and the Coulomb energy are 
comparable. The true value U = Uc of critical interaction depends on the method selected and the density of states 
for given system, but it is in the regime W < Uc ^ 2W [4]. In the limit U << 14^ we have the so-called metallic 
limit. The question interesting us here is the limit of strong correlations U » W. In the case of n = 1 this limit 
corresponds to the Mott insulating limit [2,4,5]. The question which was posed by us, I believe for the first time [6-7], 
was what happens when U » W and 1. 

The situation for n ^ 1 (it is sufficient to consider the case n < 1 only, as that with rt > 1 is related to the former 
by switching to the hole language). Namely, the Hartree-Fock approximation for U ~ 2W we have from ([4]) 




Taking n = 1 — 5, with 5 << 1, we have that Eg ~ WS/A i.e. the energy is positive and thus the metallic state is 
unstable. We show below that this is indeed the case for S < 5c 0.1, but that situation requires a refined analysis. 
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In any case, if we take Eq. ([5]) literally then for every n we can define Uc at which the two energies cancel out 
each other i.e. Uc/W = {2 — n)/n, {Uc increases as n decreases). This trivial reasoning teaches us one important 
thing: the creation of holes in the Mott insulator makes particles effectively less correlated for fixed value of U. In 
other words the condition U >> W valid for a half- filled band system gets softened when we change (diminish) 
the number of electrons in the same band. In our derivation of the t-J model [6,7] and in subsequent papers it is 
assumed that the particles are strongly correlated in the whole filling range considered. Note however that if we take 
n — 0.7 then Uc/W = 1.86, i.e. the critical value is almost doubled, so if we had band states with the value of 
U such that U/Uc — 1 for n = 1, for n = 0.7 the particles behave as if they experienced only half of the critical 
value (are substantially less correlated). This must be remembered particularly in the context of high-temperature 
superconductors described via a single-band model since then for n = 1 we have U/W ~ 2 3 [8]. 

We provide next our original version of the derivation of t-J model [6,7] and dwell on the details, which seemed to 
us not obvious whatsoever and are usually overlooked in the textbooks or review articles. 

III. T-J MODEL AS IT WAS DERIVED THEN 

I became interested in the Hubbard model around 1974 an particularly, in the problem of deriving its ferromagnetic 
state in the limit U » W . This was because on the one-hand we have a definite result of Anderson from 1959 [2] that 
the Mott insulators are antiferromagnetic and from the other that of Nagaoka [9] that the ground state of the Mott 
insulator in the limit J7 = oo is for certain lattices (e.g. fee) ferromagnetic. The obvious thing was to combine the two. 
In that manner, the question how to generalize the Anderson kinetic-exchange Hamiltonian to the metallic regime 
has arisen naturally, but the question was how to do it in detail. As the kinetic exchange is ~ /U , one should get 
it from the second-order virtual hopping processes, with the doubly occupied site configurations in the intermediate 
state. For that purpose, one can rewrite the dynamical processes contained in the hopping term as follows 

"^ia^j-y = Ac (1 ~ + nia)a.ja (1 - 71^-^ + Uj^) = (1 - ni-ff)aja- (1 - rij-ff) 

+ aj^ (1 - ni-ff)aj„nja: + a\^ni-jjaj„ (1 - rij-^) + a]^ Uia^ aj^ fija^ . (6) 

The consecutive terms represent the four restricted types of hopping processes: the first hopping from singly 
occupied site onto an empty one, the second hopping from singly occupied site and formation of a doubly occupied 
one, etc. In such manner, the full hopping term looks like 

+ ^2 ' "^'^ "i'^ ^ "^'5') + X! ' "'^'^ "^^^ ■ ^'^^ 

ija ija 

But then, there is a problem: all the processes have the same coupling constants (are coming from the same term 
in the Hamiltonian, so singling out only the second and the third term as perturbation does not sound right. After 
a few-month deliberation I will not dwell upon, I have realized that the canonical perturbation in the version by 
Bogolyubov is the right method to follow. This is because of the three following features of the method: 

1. You can perform a perturbation in an invariant (operator) form (not only on matrix elements), 

2. If iJ = Hq + Hi and Hi is a perturbation part, then in order to do the calculations Hq does not have to be 
diagonal, and you can select the perturbation term according to your wishes, 

3. It bears a strong resemblance to the quantum mechanical perturbation approach when Hamiltonian splits into 
(large) diagonal blocks and (small) off-diagonal blocks, as I learnt earlier from the Stevens lectures [10]. 

Formally, one can write the full Hamiltonian as operating in four subspaces 

Ho - ^ PoHPo + PnHPn + PoHPn + Pn H Po (8) 

ija 
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assuming (intuitively at this point) that Pq + Pn = 1, and that 

PaHPo = Y,' (1 - n,^) (1 - n,^) , (9) 

ija 

Pn H Pn = ^ ' tij Ui-ff Qjcr rij-^ + U rui , (10) 
PqHPn = ^'Ujal^il - ni-a)ajanj-^ , (11) 

ija 

PnHPq ^Y.' *y ^J'^ (1 - = (-Po ^ -P^)^ ■ (12) 

The exphcit form of these projection operators were not specified then and they need not to be. Understanding 
that took me some time. Since every term contains process ~ tij, it is convenient to define the operator 



Obviously, H = H^=i and 



H,^Ho + eHi. (13) 



Ho ^ PoHPo + PnHPn (14) 



Hi = PoHPn + PnHPo. (15) 

By introducing e we will collect the terms of the same order in e. In this manner, we will assume that PqHPq ~ t 
and PnHPjs! ~ U represent two different energy manifolds, {the Hubbard subbands) distant roughly by whereas 
PqHPm represents rare hopping process of creating a doubly occupied site from two neighboring singly occupied sites 
< z,j >. So, even though all the hopping terms contain the same hopping matrix element t^, PqHPq will provide 
contribution tij while PqHPn will provide the contribution of ~ tij /U in the first nontrivial order. 

The last nontrivial feature is to select the unitary transformation matrix, with the help of which we remove part the 
hopping processes in the first order and replace them by virtual processes exemplifying physical processes in higher 
order. For that purpose, we have also proposed the following canonical transformation of the form 

H, = e-^'^He^'^ ^ Hq + e (iJi +t[Ho,S]) - i {2i[Hi,S] - [[H^ , S] , S]) (16) 

with S = S'<. 

The linear term ~ e is absent when 

Hi + i[Hq, S] = (17) 

Under this condition 

He ^ Ho + ^ie^[Hi, S] (18) 

Eq. (|17p has to be solved for 5. The main problem is that Hq is not diagonal. So, it is more difficult to solve it 
the way it was treated originally by Bogolyubov (cf. also the treatment of electron-phonon case by Kittel). Thus, in 
order to proceed one tries to solve (|17p projected onto the subspaces. Explicitly, we can write 



PoHiPn + iPo[Ho, S]Pn = 0, 



(19) 
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Po {Ho S - SHo)Po = Pn {Ho S - S Ho) Pn = . (20) 

From the first equation we see that while S must have nontrivial form as it contains Hi , the next two equations are 
of trivial character. This is because one can rewrite them in the form 

(Po Ho Po) {Po SPo)- {PoS Po) Po {Ho Pq) - , (21) 



(Pat Ho Pn) {Pn S Pn) - {Pn S Pn) {Pn Ho Pn) = . (22) 

In the original version we have taken a particular solution Po,n S Po^n — Jo,nPo,n- Equally good would be to replace 
r.h.s. with an arbitrary function /(Pq^at) which has a Taylor expansion as /(Pq^at) — ao.N^ + ^o,nPo,n because of 
the property P^ j^ ^ Po,n- 

In effect, we are left with Eq. (fT7|) . which projected takes one of the two possible forms 

PoSPn = i-i Po Hi Pn + {Po H Po) {Po S Pn)] {Pn H Pn)'' . (23) 

Obviously this relation is well defined if the operator {PnHPn)~^ exists, but we leave such problems to mathematical 
physicists. What is important for us instead is, that since the transformation matrix appears on both l.h.s. and 
r.h.s. of we can try to solve it by iterating the solution. In the zeroth order we assumed that on r.h.s. of (|^ 

{PoS'-"^Pn) = and then 

Po ^(1) Pn - -^ Po Hi Pn I Pn H Pn . (24) 

The solution up to infinite order takes the form 

Po ^ Pjv = -I Po Hi Pn {Pn H Pn - Po H Po)-^ . (25) 

The simplest approximation of the denominator is to replace it by an average < PnHPn — PqHPo >~ U , where U 
is the energy difference between the centers of gravity of the Hubbard subbands. In effect, after substituting ([25|l to 
(fT8|) as well as its Hermitian conjugate PnSPo, we obtain the effective Hamiltonian H = H^^i up to P/U as 



H ^ PoHPo + PnHPn, (26) 

with 

PoHPo = PoHPo - PoHPnHPo/U, (27) 

and 

PnHPn - PnHPn + PnHPoHPn. (28) 

The first part of ((27|l describes the dynamics of electrons in the lower Hubbard subband for n 7^ 1, whereas (f28| will 
represent the same for 1 < n 7^ 2. In the following discussion we limit ourselves to the situation with n =/= 1. 

A. Meaning of the t-J model 

Substituting (fTTj) and (fT2|) to (p7)) and (|28|) . respectively we obtain exphcitly [11] 



PqH Po = ^ Uj (1 - Ui^) Qja (1 - Uja) 
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■ni-g)n^„f (1 



+ 



E 



(1 



(1 - riior) S ■ flfcs^ (1 - rifco-) 



(29) 



and 



N 



^i(T ^"io" '^jcr ^JCT 



2^2. 

^ «!t "a"ji + . . . . (30) 



The first term in (|29p represents electron hopping from single-occupied site onto an empty site i, the second term 
represents full kinetic-exchange part in which, in general, the spin operators {Si} representing itinerant electrons are 
expressed in terms of fermionic operators: Si = {S^ , S~ , Sf) = {a\^ai]^, a^^ , (n^i — 12^)/ 2). The third term denotes 
hopping between three sites without and with spin flip, respectively, in the middle (j) site. One can easily see that in 
the limit of Mott-Hubbard insulator when the number of particles is conserved and equal to unity at each site (i.e. 

+ nil ^ ^'-'^ only < ti^i -|- >= l!), then (j29p reduces properly to the Anderson (antiferromagnctic) kinetic 
exchange Hamiltonian 

PoHPo = V ^ fs, • S, - i V (31) 



^ U V ' ' 4 

In our original work [7,6] we have already noticed that the dynamics requires introduction of projected fermionic 
creation, annihilation, and particle-number operators 

bl ^ al (1 - n,^) , (32) 
bia = flifT (1 - rii-^) , (33) 

i^ia = bi„ = riio- (1 - riia) , (34) 
which have nonfermionic anticommutation relations 

= [(1 - - S-'^il - S^^,)] (35) 



{h.,b,^,} = {bl,bl,} =0 (36) 

The property (j35p poses a problem as there is no reference occupation-number representation. Additionally, the 
hopping and the exchange parts do not commute with each other, so the motion of single electrons and the spin 
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interactions (e.g. spin flip processes) are entangled with each other and result in a strongly correlated (non-Fermi 
liquid) metaUic state. What is equally important, the kinetic energy (~ tij) and the exchange {tfj/U) can become 
comparable for a filling close to the half filling and we can have a transition to the localized state induced by the 
exchange interaction (see next Section) or to the phase separation into antiferromagnetic insulating islands " fioating" 
in the ferromagnetic sea of (1 — n) holes [9]. This second state is rather unstable as the Coulomb repulsion introduced 
by the holes i.e. by the positively charged ions of the background (usually neglected in the calculations!) will 
destabilize the phase-separated configuration. 

Few comments are relevant at this point. First, the operators were introduced by Hubbard [12] under the 

name of atomic representation. The original Hubbard notation is useful for multi-orbital models [13]. Nonetheless, it 
is used sometimes also for the present model though, in our view, represents an unnecessary formality. Second, it is 
sometimes important to include the intersite Coulomb interaction, i.e. the term 



In that situation [9], the kinetic- exchange integral takes the form Jy = 2tj^/{U — Kij) and an additional term 
{1/2)^^-' Kij ViVj appears in (|29p . Finally, the essentially the same formalism has been repeated later by other 
authors [14]. It is unfortunate for us, those references are often quoted as those, from which the t-J model (i.e. (|29p . 
without 3-site terms present) originated. We also comment on the effective t-J model later on when we introduce the 
two-orbital model. 



B. Localization of holes in doped Mott insulator 



Intimately connected with the t-J model is the problem of magnetic polarons [15] and the hole localization induced 
by the kinetic exchange for the small number of correlated holes, 5 — l — n<^ 1. One can address the hole localization 
by starting from the metallic side. The ground state energy of ([29]) can be estimated as follows 

Eg^ {H) ^t J2 («L(1 - "-)aja(l - n,^)) + J ^ (^^ ' " l"^"^) ' (3^) 

<ij>a <ij> 

where we have included the terms involving only z nearest neighbors < i,j >. In the spirit of Gutzwiller approach, 
we can renormalize the hopping term by the band narrowing factor $(A), with A = {viVj/A — Sj • Sj)/n'^ and thus 
have 



= zt^X) J2 {ala,,)^~ Jzn'X (39) 

where (. . .)o denotes the average for an uncorrelated state. Next, expand $(A) — go + giX + 52A^, while {al^aja)o = 
7icr(l — n) expresses the hopping probability in the reference state without antiferromagnetic correlations (A = 1/4). 
For that value of A we have go + 5i/4 + (72/I6 = 1. In the opposite limit (A = 1), we have a complete frozen 
antiferromagnetic (Neel) state, without any hopping, i.e. go + 5i + 52 = 1- Furthermore, writing the energy in the 
form 



= -\t\n{l - n){go + gi X + g^X^) - J X , (40) 

z Jy 

and minimizing it with respect to A, we obtain the expression 



4g2A = -^- r+ffi. (41) 

|t| (1 - n) 

This equation must also fulfilled for J = 0, when there are no antiferromagnetic correlations (i.e. when A — 1/4). 
This means that from the above three conditions we can determine explicitly the coefficients of the expansion of *I'(A), 
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which have the values: go = —1/9, gi = 8/9, and 92 = 16/9. Hence, the optimal value of the variational parameter A 
is 

1 1+9J 

This means that the full antifcrromagnctic insulating state (A = 1) is achieved for 



S-Sc^[l + fj) , (43) 

and is of the order Sc — 0.1 for \t\/J = 3, a reasonable value for La2_2;Srj;Cu04 system. The ground state energy is 
then 

Eg , , . , 1 2 9 J^n^ 

-| = - t n{l - n) ~ -Jn' - - — ^ . 44 

z Jy 4 64 \t\ (1 — n) 

This reasoning represents essentially the simplest mean-field-like description of the transition from the strongly 
correlated paramagnetic metal to an antiferromagnetic insulator with frozen holes at S ^ 6c- 

IV. T-J MODEL RENEWAL: KINETIC EXCHANGE AS THE SOURCE OF REAL SPACE PAIRING 

Until 1987, the t-J model was regarded as a generic model for explaining the antiferromagnetism of Mott insulators 
such as the oxides NiO or CoO, and relatively rarely considered for n 7^ 1. Obviously, the real oxides have both 3c? 
and 2p orbital states as the valence-band states, but it was argued strongly (and still is) [2,15,16], that the p — d 
hybridization can be incorporated into an effective d-band Hubbard model, although there were some exceptions to 
that idea [17]. The essential new idea, about which I have learnt from the preprint of Ruckenstein et al. [18] was to 
associate the kinetic exchange part with real space pairing, but there were quite few papers appearing at about the 
same time [19]. Soon, the slave-boson and slave- fermion mean-field theories have been formulated followed by a flood 
of papers, also about the stability of d-wave superconducting solution. This is a fascinating story for itself, but I shall 
leave this to a separate occasion. 

A. t-J Hamiltonian as pairing Hamiltonian 

I would like to raise only one aspect of this story namely, the proper formal expression of the pairing part [11]. As 
we have seen, we have to express the whole Hamiltonian (1291) in terms of projected fermion operators hic^ and b\^. 
First, one notices that the spin operator has the same form in both fermion {0^0-, a\^} and in the projected {6^0-, 
representations. Explicitly, 

= (^a\^ a,x , o^t , (n.-r - ni|)/2^ = {h\^ , , (»^iT " '^a)/2) ■ (45) 

Therefore, the spin-singlet real-space paring operators should be selected in the form (i 7^ j): 



Bl 



iK^li-^Mr)' (46) 



V2 



Bij = {bi]bji - hibj^) ^ --^ (aiTOji - aiiaj]) (1 - n^^^) [l - n^i) . (47) 



In this representation, we have the following identity 
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4 Brj ^-[^^■^3-\Y1 ■ (48) 

One should note at this point that the representation of (|48| through the operators bI- and Bij is completely 
equivalent to that through the spins and projected particle-number operators. In other words, the spin ordering and 
the nearest-neighbor pairing should be treated on the same footing. 

Substituting ([^5]) to PqHPq we obtain the effective t-J Hamiltonian in the following closed form when the 3-site 
terms are also included [11] 

H ^ PoHPo - J2 - E 2 tjk/U 4 B,k ■ (49) 

ija ijk 

Obviously, it is assumed, that ty ^ only for i ^ j. Furthermore, the pairing part vanishes by itself if i = j as 
bI = 0. This means that the pairing order parameter will have a nontrivial dependence in both real and reciprocal 
(k) spaces. 

Transforming the Hamiltonian (|49p to k space, we obtain pairing Hamiltonian, which takes the following form when 
there are phase differences between the neighbors from the left and the right sides. 

PoHPo = E 6L ^k. - ^ 5] 7k 7k' 4'-k ^-k'k , (50) 

ko- kk' 



where 
and 



4'_v ^ ^ {bl &-ki - bl^ 6_kT) (51) 



5L = 47E^"^'''''«l(l-"^-) (52) 



From the form ((52)) we see that also the operators b]^^ (and feko-) have a nontrivial many-particle character. The 
only simple property here is the separability of the pairing potential, i.e. Vkk' = ^(2t^/?7)7k7k' = VkVk', where 7k 
depends on the solution symmetry (for example, for an extended s-wave solution for the square lattice 7k = ^k/t] see 
also below). 

The solution of t-J model in the forms or ([50)1 is the subject of an intense discussion, but this absolutely 
fundamental aspect of the model will not be touched upon here. We hope to return to this aspect of the exchange 
mediated superconductivity elsewhere. 

B. Three remarks on the pairing Hamiltonian with 3-site terms 

First remark concerns the formal form of 7k. Namely, when taking the Fourier transform of the part tijtjk in (I49p 
and assuming that for given neighbor (e.g. in the square lattice) the operator Bij for j to the right from i acquires 
the phase (px with respect to that to the left. In effect, the part 7k of the pairing potential can be expressed as 

— gifcxO-l-Vx _|_ g-ifcxfl _j_ ^ikyU + ipy _|_ ^-ikya ^gg-j 

or equivalently 

7k = 2 cos (^kx + ^) cos (^^j -I- 2i sin (k^ + ^) sin 
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+ 2 cos (ky + ^) COS (^) + 2 i sin (fc, + ^) sin (^) (54) 

This form represents the most general k dependence of the pairing when the 3-site terms are included. 

Second remark concerns the equivalence ([48]) of spin degrees of freedom and the real-space pairing operators (note 
however, that r.h.s is nonzero for i = j). Therefore, when decoupling r.h.s of (|^^ we should take into account a 
possibility of have nonzero both averages. In the Hartree-Fock-like approximation this decoupling for k — i should 
take the form 

4 B,, ^ (4> B,, + {B,,) 4 ~ (4) {B,,) - [5f + 5f - (5f> (S^) ] + . . . (55) 

So, in general, the coexistence of antiferromagnetism (or spin density wave) and the paired state is possible and 
should be studied carefully. 

Last remark concerns mapping of the projected Hamiltonian onto an effective fermionic Hamiltonian, as proposed 
[11,20] on early stage of the theory development. Namely, introducing the concept of band narrowing $ = <i>(A) one 
can have that 

H - (1 - $) 5] {2tllU) Bj^ B,, 

ija ij 

- $ (1 - $) ^ {2U,tk,/U)Bl Bk, , (56) 

ijk 

where < $ < 1 is the band narrowing factor. Macroscopically $ can be regarded as a degree of itineracy of 
electrons. Likewise, (1 — $) is the "fraction" of localized electrons. In ([55]) the two-site and three-site terms have 
different renormalization factors. Hence, is tempting to take, as earlier, {B\-Bij) as the variational parameter A. 
Then, as before $ = $(A). It is also tempting then to associate the pairing only with the 3-site part, 

but as we are aware of, this possibility has not been studied explicitly as yet. In such model the limit of the 
antiferromagnetic insulator is reached naturally for (5 — s- and the pseudogap appearance is associated with short- 
range antiferromagnetic correlations. 

V. KONDO INTERACTION AND REAL-SPACE HYBRID PAIRING AS AN EXTENSION OF THE T-J 

MODEL CONCEPTS 

A. Deep and shallow impurity cases 

Shortly after deriving the effective Hamiltonian and determining the principal features of antiferromagnetic 
phase I have turned my attention to the Wolff model of magnetic impurity [21] and the corresponding generalization 
of the canonical transformation [7,6], as it has been proposed and discussed by us earlier [22]. The Wolff model differs 
from the original Anderson-impurity model [23] only slightly, as it treats both the impurity electron state and the 
band states of a metal, the impurity is immersed in, in the tight binding approximation. Explicitly, the Hamiltonian 
of the system is proposed in the form 

H = ^"iy aLojo + e<i("0T + '^oi) + U not noj 

ija 

+ (1 + 7) XI (°0o + a]^ ao,,) . (57) 

JO" 

Here we assume that the impurity is at the 0th site, with as the energy for the one-electron impurity state. The 
intraatomic Coulomb energy for the two electron impurity state is U (and zero for the metallic host sites J 0). The 
hopping strength between the impurity and the host is specified by 7. The double primed summation in (j57p excludes 
the terms with either i = or j = 0; is the impurity energy-level position with respect to that for band states, 
which is taken as e,, = 0. 
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We were interested in the strong correlation limit when |(1 + 7)toj"| is much smaller than U. What was (and still 
is) new in our original approach is, that we have considered separately the shallow- and deep-impurity limits. This 
turned out to have a direct relation respectively to the Kondo and to the mixed- valence (or heavy-fermion) limits of 
correlated electrons. Also, it has relevance to the concept oi p — d pairing in the high-temperature-superconductivity 
era [24]. As these reminiscences are probably too long-winded already, I provide only the essence of the argument 
leading to the pairing. But first, summarize the results for the impurity situation. 

In the interesting us case of a shallow impurity wc have the impurity (correlated) state immersed in the Fermi sea 
of the metal. In that situation singly occupied impurity state is relatively close to the metal Fermi energy ep but its 
doubly occupied state is high in energy and therefore, can be realized via virtual excitations (hence, transformed out 
again). In effect, the effective Hamiltonian in the second order takes rather complicated form 

H = PqH Po = Y" tij al^Uja- + (1-^7)^ toj (1 - now) (al^aja- + o],^ ao^) 

ijrr jrr 

+ Q (^OT + ^oi) + ^^^^ J^' J2 ^oj ( So • S, - i ^ i^oa ) • (58) 

j \ crcr' / 

This effective interaction contains both the antiferromagnetic impurity-host kinetic exchange interaction (the last 

term, which can be called the non-local Kondo interaction) and the residual hopping between the impurity and the 
metallic host. This form, at first look, seems to be of limited use, as there is still part of the hopping left between 
the host and the impurity. As we show in the following, this type of situation is specific for the heavy fermions and 
high-Tc systems, since the valence {p) states play an active role in the electric conductivity. 

For the sake of completeness, we provide also the effective Hamiltonian in the deep-impurity case, i.e. when the ea 
level is deep below the Fermi level. It is 



2(1 + 7)^^ 
ed {U + ea) 




^ ^ Uo tjO Sq al- Qja- 



(59) 



This Hamiltonian has a direct correspondence to that derived by Schricffer and Wolff [25] for the Anderson-impurity 
model. The last term represents the Kondo interaction, with three-site spin-flip term included (here we have explicitly 
assumed, that in the present limit not + "-oi = i-^- we have localized moment at the impurity (0th) site). 

B. Hybrid pairing in two-orbital system 

After the t-J model had been interpreted in the categories of electron pairing, a natural question arose if one can 
do the same for the hybridized (two-orbital) system such as the Anderson-lattice model. This idea has been put into 
writing subsequently [24] by starting from the Hamiltonian 

H = 'tmncl^aCna + Cf^^i'^ + ^Y 

mna ia i 

+ Y + ^^rna "-i^) ' (60) 

im 

where the {i,j) label correlated atomic (o = d or /) states, {m,n) label delocalized c states, and Ni„ = a\^aia- Vim 
represents hybridization integral. As in the case of shallow impurity, we assume that \Vim\ -C U, but |e/| ~ \Vim\, so 
not the whole hybridization term can be transformed out. Now, decompose the hybridization term as before, i.e. 
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al^Cnia + C^raa^'i'y = (1 ^ ^ ia) {a\a C,na + c\^„ ai„) + Ni„ {a\^ Craa + C^raa ^'i'y) ■ (61) 

The first term on r.h.s. represents the interstate (a — c) hopping processes which do not involve U (double occupancies 
of atomic states), while the second does involve U and hence leads to higher-order mixing processes. The basic idea 
introduced at this point [11] is to canonically transform out the second term only and to replace it by an effective 
interaction incorporating higher-order virtual processes. Leaving the details [24], the effective Hamiltonian with 
real-space pairing has the form 



^ — ' U + Cf 

imn 

where the hybrid-pairing operators are defined as 

Bin ^ ^ («It ^li - «a ^U) (1 - • (63) 
This model can be transformed again to the effective Fermi-liquid form and the proper form is: 



2 



U + ^f 

where the tilted quantities are renormalized by correlation, e.g. 



E^kl^'l^-k^q.-q^ (64) 



Vi. = Vk - V e**"^^ Vjrn , (65) 



j(m) 



and = (1 — «/)/(! — "•/o-), with n/^ = (alg-aitr). Also, the factor (1 — N^-ff) is absent in (cf. Eq. ([63]) ). This 
Hamiltonian can be solved in BCS approximation, but the results will not be reproduced here. Likewise, we will omit 
here the application of (|62|) to the description of high-Tc superconductivity with hybrid p ~ d pairing, as it would 
most probably require the inclusion of the 4th order (d — d or / — /) interactions. 



C. Kinetic exchange in orbitally degenerate systems and possibility of spin-triplet pairing 

This topic grew into an independent discipline of its own after publication of the papers by Kugel and Khomskii [26] , 
Lacroix and Cryot [27], and Inagaki [28]. Some aspects of the topic are reviewed in this issue by Oles [29]. We have 
also derived the kinetic- exchange Hamiltonian for a partial filing of the band [30] , as well as have applied it to explain 
ferromagnetism of C0S2 [31]. In general, inclusion of the orbital degeneracy allows for a natural explanation of the 
appearance of ferromagnetic Mott-Huhbard insulators in conjunction with the orbital ordering of antifcrromagnctic 
type in systems such as K2CUF4. 

In connection with this, a question has arisen, whether the ferromagnetic kinetic exchange can produce spin-triplet 
pairing. This question is particularly important because it has been shown earlier, that the Hund's rule coupling 
can lead not only to ferromagnetism, but also to the superconducting pairing [32]. Such an effective model has been 
proposed and analyzed in detail in the strong-correlation limit very recently [33] . These last results will be submitted 
for a publication shortly. 
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VI. CONCLUDING REMARKS 



In this overview we have concentrated on the idea of kinetic exchange as it was derived 30 years ago and its 
subsequent application to the systems with real space pairing mediated by this exchange interaction. By no means 
it is a complete survey. For example, we have ignored all the subsequent formal development of the model (see e.g. 
[34]). Also, we have disregarded the effect of electron-lattice coupling on the effective t-J model with pairing [35]. 
Nonetheless, what I hope I have sketched here is the analytical structure of the t-J model in various strongly correlated 
systems, in which magnetism and superconductivity seem to have a common origin, although it has not as yet been 
proved conclusively, that the kinetic exchange is the origin of both of them. Future will show. 
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